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Abstract. An optimization problem is formulated motivated by the de- 
sire to remove temperature peaks, i.e., “hot spots”, along the bounding 
surfaces of containers of fluid flows. The heat equation of the solid con- 
tainer is coupled to the energy equations for the fluid. Heat sources can 
be located in the solid body, the fluid, or both. Control is effected by ad- 
justments to the temperature of the fluid at the inflow boundary. Both 
mathematical analyses and computational experiments are given. 


1. Introduction 

We suppose that the regular bounded domain Q, in R 1 2 is made up of two sub- 
domains fli and separated by an interface IV with the result that Q. = fliUftjU 
T w (see Figure 1). The solid material occupies a sub-domain having a boundary 
r x U r 2 u r 3 U IV and the fluid flow occupies a domain Cl 2 having a boundary 
r 4 U r, U r„ U r 4 . We have an inflow boundary IV an outflow boundary IV and 
a solid wall IV The geometry of all these boundary segments is prescribed, as are 
the inflow velocity u e and temperature T e . At the outflow, one can impose one’s 
favorite outflow boundary conditions. On the walls, we have the no slip boundary 
conditions for the velocity. Control is to be effected through heating and cooling 
along the boundary IV 

The temperature is specified along the boundary T e and the heat-flux is specified 
along the boundary I\ U IV U T 3 U I? 4 U IV We assume that the flow is incompressible 
and convection driven so that buoyancy effects can be neglected, and thus tempera- 
ture effects on the mechanical properties of the flow, i.e., the velocity and pressure, 
are negligible. We are interested in controls such that we get a desired temperature 
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r 4 

Figure 1. The Domain 12 


along T w or a portion r„ C r w , and thus we assume that the flow is stationary. 
Other combinations of control and controlled surfaces are also possible. 

As a result of our assumptions about the flow, the state variables, i.e., the velocity 
u, pressure p, temperature T, and control g are required to satisfy the Navier-Stokes 
equations 

(1.1) - uAu + (u • V)u + Vp = f in 12 2 ? 
the incompressibility constraint 

(1.2) div u = 0 in 
and, for simplicity, the boundary condition 

(1.3) 

(1.4) 

(1.5) 

and the energy equations 

(1.6) — Ki&T = Qi infix. 


u = 0 on Tu, U r 4 , 
du 

^ = 0 onr„, 


(1.7) - k 2 AT + (u • V)T = Q 2 + 2/r(Vu + Vu T ) : (Vu+ Vu T ) in fi 2 , 
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with the boundary conditions 

(1.8) T = g on T e , 

, . 0T 

(1.9) ~dn =0 on Ti u r 2 u r 3 u r 4 u r 0 . 

The data functions f, h, Q x and Q 2 are assumed to be known. The constant u is the 
kinematic viscosity coefficient of the fluid, and the constants #c lt k 2 and p. depend 
on the thermal conductivity coefficient, density, specific heat at constant volume, 
and viscosity coefficient of the fluid; see [17] for details. 

Note that as a result of our assumptions about the flow, the mechanical equations 
(1.1)-(1.5) uncouple from the thermal equations (1.6)-(1.9). Indeed, (1.6)-(1.9) may 
be solved for u and p without regard of the temperature T. Thus, in the present 
context, the velocity field u, which is determined by solving (1.1)-(1.5), merely acts 
as a coefficient function and in the source term in (1.7). 

We now define the optimal control problem to obtain our objective. For example, 
given a velocity field u, we would seek a temperature field T and a control g such 
that the functional 

(1.10) ir-iif (Ijf+iv.jivr 

is minimized subject to ( 1.6)-( 1 .9), where V, denotes the surface gradient operator 
and T d is some desired temperature distribution, e.g., something close to the average 
temperature along r„ for the uncontrolled system. The non-negative parameters 7 
and 6 can be used to change the relative importance of the two terms appearing 
in the definition of J as well as to act as penalty parameter. Incidentally, the 
appearance of the control g in the J is necessary because we are not imposing any 
a priori limits on the size of this control. 

Under the realistic assumption that u • n = 0 on U T 4 and u-n > Oon f 0 , 
in this paper we prove the existence and uniqueness of optimal solutions and drive 
an optimality system, i.e., a set of equations from which the optimal control and 
state may be determined. Also, a finite element method is used to compute an 
approximate solution of the optimality system. We have also developed an iterative 
algorithm to compute the approximate solution. 

We close this section by introducing some of the notation used in subsequent 
sections. Throughout, C will denote a positive constant whose meaning and value 
changes with context. Also, H’{V), s € R, denotes the standard Sobolev space of 
order s with respect to the set X>, where V is either the domain ft , or its boundary 
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r, or part of that boundary. Of course, H\V) = L 2 {V). Dual spaces will be 
denoted by (•)*. 

Of particular interest will be the space 

(1.11) = {S6 ^ G £ 2 (ft) for k = 1,2} 

and the subspace 

(1.12) = {5 6 H\a) : S = 0 onTj. 

For functions defined on r c , we will use the subspace 


(1.13) W(T e ) ={g£ H'iTc) \ g = 0 at f e n f,}. 


Norms of functions belonging to H’{ T) and H’(T e ) are denoted by [| ■ Ha.n? 

|| • ||, r and || • || , respectively. Of particular interest are the I 2 (Q)-norm || • || 0 ,n> 
the semi- norm 


(1.14) 


\T\l « 



2 

o 


and norm 

(1.15) imi 2 , n = |T| 2 , n + ||r|| 2 )n 

defined for functions belonging to H\Q). Also, we are interested in the semi-norm 
| • |i.r e , defined by 

(1.16) lsl?.r c = jf IV.pl 2 ^ 
and norm 

(1.17) Il5lli,r c = l5li,r e + I l^llo,r c 

defined for functions belonging to H l {T c ) and W(r e ). 

We define, for (T5) 6 £*(11), 

(1.18) (T ’ 5)n= Ijsdn 


and, for ( pq ) € L l (T), 

(1.19) (p,q)r = j^pqdT. 

Thus, the inner product in £ 2 (ft) is denoted by (•, -)n, that in L 2 (T) by (-, -)r- The 
notation of (1.18>(1.19) will also be employed to denote pairing between Sobolev 
spaces and their duals. 
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We will use the bilinear forms, for i = 1,2, 

(1.20) ai(T,S) = Ki[ VT-VSdft VT,S € H l (Q) 

Jtii 

with the result that 

^ r 

(1.21) a(T, S) = 5> / VT-VS<m = ai (T, S ) + a 2 (T, S), 

.=i 

and the trilinear form 

(1.22) c(u,r,5)n a = [ (u-VT)SdQ Vu G and VT,S e H\n 2 ). 

These forms are continuous in the sense that there exist constants c, > 0 and c c > 0 
such that, for i = 1,2 


(1.23) \<H(T,S)\ < c I ||T|| 1 , n ,||5|| 1 , n . VT,S G F l (fl), 

(1.24) |c(u,T,S)n 3 | < c # |N| lA ||r|| ll o t || 5 || ll0l 

Vu G H 1 (fl 2 ) and VT, S G 


and thus 


Kr, 5)| < mt, 5)|+ Mr, s)\ 

(1.25) < c 1 ||T|| x ,„ 1 ||5|| 1 . ni + c 2 ||r|| 1 , n J|5|| 1 .n a 

< (ci + c 2 )||r|| li n||S||i > n VT, 5 G H\n). 

Moreover, we have the coercivity property, for i = 1,2, there exist constants C,- > 0 
such that 

(1.26) o,(T, T) > C<||r||J n . VT G ‘(ft) 

and thus 

(1.27) a{T,T) = a t (T,T) + a ,(T,T) > £C,||T||? i0l > min(C 1 ,C I )||r||f,„. 
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2. The Optimization Problem and the Existence of Optimal Solutions 

We begin by giving a precise statement of the optimization problem we consider. 
We will assume the domain f2 is in R 2 and consists of two subdomains fij and 
such that U Q 2 U IV Let g € W(T e ) denote the boundary control and 

let T G H l (Q) denote the state, i.e., the temperature field. The state and control 
variables are constrained to satisfy the system eqrefeq:asstl-(1.9) , which we recast 
into the following weak form: Find ( T,t ) G 5 rl (fi) x jy -1 * 2 (r e ) such that 

(2.1) a{T,S) + c{n,T,S) n3 -(t,S)r c = (Q,S) n VS G H 1 ^) 
and 

(2.2) (T, R)r c - ( 9 , R)r c = 0 ViZ G E~^ 2 (T e ), 
where we have introduced the simplifying notation 

n _ f Q 1 in 

^ \ Q 2 + 2/x(Vu + Vu r ) : (Vu + Vu T ) in ft 2 . 

One may show that, in the distributional sense, 

(2.3) t = k 2 VT • n| r< _. 

In (2.1)-(2.2), we introduced the Lagrange multiplier t to enforce the boundary 
condition. This will be very useful in the proof of error estimates for finite element 
approximations. 

First, we show that for each possible control function g, there is a unique corre- 
sponding state function ( T,t ). 

Lemma 2.1. For every g G W(r c ), there exists a unique (T,t) G 
such that (2.1)-(2.2) are satisfied. Moreover, there exists a constant C = C(ti) such 
that 

(2.4) ||T||i,n + ||<il-i/ 2 , r e ^ C(||Q||-i,n + IMIi,r c )- 
Proof. For given g G W(r c ), (2.1)-(2.2) is equivalent to 

(2.5) a(r,S) + c(u,r,S)n, = (Q,S)n VS G 

T = g on T e 


( 2 . 6 ) 

and 

(2.7) 


t = k 2 VT ■ n| Fc - 
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By the virtue of the trace theorem, let T in H 1 ^) satisfy f = g on T e and examine 
the following problem: Find T in H l (Si) such that 

T-T £ 

(2.8) a(T — T,S) + c(u,T—T,S)a 3 

= (Q,S) n -a(T,S)-c(u,f,S) n a V5 € ffi(ft). 

Let T = T — T £ From the assumption of u, t.e., u • n = 0 on r u U r 4 and 

u • n < 0 on r 0 , we have that 

c(u,f,T)n I = J / (u • n) T 2 dT 

(2.9) ^ • /d ° 3 

= -^(u-n)T 2 dr>0. 

Thus, we have 

(2.10) «(T, f ) + c(u, r, r)n, > min(C 1 , C 2 )||f || 2 _ n . 

Therefore, by the Lax-Milgram theorem there is a unique T £ Hp(Q), t.e., T = 
T + T £ 5’ 1 (Q) and the estimate 

(2.11) imu < C(||g||. 1 ,n+ |Mli/ 2 .rJ < C'(||g||_ 1 , n + ||*|| 1>r .) 

holds. From the trace theorem and the theory of partial differential equations (see 
[4]), we have 

(2.12) IWI-i/2, r, < CdlTIU.n, + IIQnll-i.n,) < C(||T|| 1>n + \\Q\\_ h n), 

where Q 2 = Q? + 2/x(Vu + Vu T ) : (Vu + Vu T ). Then, (2.4) follows from above two 
estimates. □ 


The admissibility set U a< i is defined by 

H a d = {( T,g ) E x W(T C ) : J(T,g) < oo,and there exists a 

t E H~ l ^ 2 (T c ) and (2.1)-(2.2) is satisfied}. 

Then, ( T, #) E is called an optimal solution if there exists e > 0 such that 

(2-14) 

J(T,g) < J(T,g ) V(T,p) 6 U ad satisfying ||T - f||i + ||<? -5||i,r e < e. 
We now show the existence and uniqueness of optimal solutions. 

Theorem 2.2. There exists a unique optimal solution ( T,g ) £ U ad . 
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Proof. We first claim that U ad is not empty. Let g = 0 and then let (T,g) G 
H'iQ.) x W(r c ) be a solution of (2.1)-(2.2); note that with g = 0, (2.1)-(2.2) is 
equivalent to 

a(f , S) + (u, T, S) n , = (Q, S) a VS G #£(«), 

T = 0 on T c , 


and 

(2.15) t = k 2 VT • n| re . 

By Lemma 2.1, (T,f) exists and (f,0) G ZV od . Now let {T (n) ,ff (n) } be a sequence in 
U a d such that 

(2.16) a(T<">, S) + (u, m S) n , - (t (n) , S) r c = (<9, S)n VS G 

(2.17) (T<">, *)r« - (5 (n) , * )r e = 0 ViZ G #" 1/2 (r c ), 
and 

lim J(T<"W n >) = inf J(r, 5 ), 

(T,^)€W a<< 

for some t (n) G #~ 1/2 (r c ). Then, using (1.10) and (2.13), we have that {||ff (n) ||i,r c } 
is uniformly bounded which in turn yields that {||T^ n ^||i} and { 1 I— i/ 2 ,r c } are 
u nif ormly bounded. We may then extract subsequences such that 

g {n) -9 in W(T e ) 

T (n) — fin H l {Cl) 
t (n) - t in ff~ 1/2 (r e ) 

T (n) - f in Z 2 (ft) 

T (n)|^ ±\^ in L 2(J g ) 

for some (f , g) G S’ 1 (ft) X W(r c ). The last two convergence results above follow 
from the compact imbeddings if l (ft) C L 2 (ft) and H 1 ^ 7 ( T,,) C X 2 (F <7 ). We may 
then easily pass to the limit in (2.16)-(2.17) to determine that ( T,g,t) satisfies (2.1)- 
(2.2). Now, by the weak lower semicontinuity of J(-, •), we conclude that (T,g) is 
an optimal solution, i.e., 
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Thus, we have shown that an optimal solution belonging to U ad exists. Finally, the 
uniqueness of the optimal solution follows from the convexity of the functional and 
the linearity of the constraint equations. □ 


Let 


2 o 2 a 2 



be a differential operator of the second order in divergence form on an open set Q. 
of R 2 . We introduce the bilinear differentia] operator associated with L 


(2.18) 


Now, setting 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 


T (T c\- ST dT dS , V' dT, 

L(T, S ) - £ a i3 + E c « dx . S - 

» J = 1 * ^ *2=1 * 


Jk 1 

on f2] 

a ll — a 22 — \ 

on f2 : 

[k 2 

a 12 — a 2 i — ■ 0 

on f), 


fo on fli, 
\ui on fi 2 , 


where u = (u 1? u 2 ), we have the following theorems. 


Theorem 2.3. LetT be the solution of (2.5)-(2,6) and letT x andT 2 be the restric- 
tions of T to fti and tt 2 , then T x and T 2 are solutions of a transmission problem 

(2.22) a x (T u S x ) = (Q U S X ) ni V5 X e JT 1 ^), 


( 2 . 23 ) a 2 (T 2 , *$ 2 ) + c(u, r 2 , 5 2 ) = (C? 2 i ^ 2)03 V 5 2 6 #i>(ft 2 )> 

(2.24) T 2 = g on T c , 

(2.25) T x = r 2 on r w , 

fW f) r P 

(2.26) = 0 ° nr, “- 
Proof. For the proof, see [12]. □ 
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Theorem 2.4. If we suppose the restrictions of Q in (2.5)-(2.6) to fij and Q, 2 are 
U r») and C°°(Q 2 U r„), respectively, and g G C°° , then every solution 
in ,ff*(fl) of (2.5)-(2.6) has its restrictions T\, T 2 to fii, 1^2 C (I^i G Tu,), 
ur„), respectively, and therefore in particular the solutions of the transmis- 
sion problem (2.22)-(2.26). 

Proof. For the proof, see ([11], Proposition 9, p592 ). □ 

We define the space H’(Q) for s > 1 by 

(2.27) H’(V) = {T G J7 1 (fi)|||r||*. (n) < 

where 

(2.28) ll^’lllr > (n) = II^IUtho) + 53 

i=i 

Theorem 2.5. Let T € JT^ft) be the solution of the problem (2.1)-(2.2). Then, we 
have 

(2.29) ll^llA^n) + ll*llir*/»(r») 5 C||<9IU°(n)- 
Proof. For the proof, see ([5], Theorem 8.5.1 ). □ 


3. First-Order Necessary Conditions for the Optimal Solution and 

an Optimal System 


We now proceed to derive the first-order optimality conditions associated with 
problem (2.14). The optimal control problem (2.14) is equivalent to the following 
minimization problem: Find g 6 W(r c ) such that K.{g) := J(T{g),g) is minimized 
where T{g) G is defined as solution of (2.1>(2.2). By studying the Gateaux 

derivative of the functional K.(g), we can obtain the first-order necessary conditions 
for the opt imal solution (T, g) in a straightforward manner. Let g be a solution of 
the minimization problem min, € »y(r e ) £($), then for ever y 2 6 W(T C ) we have 


(3.1) 

VA 

due to the definition of 

g. In 

(3.2) 

VA > 

and, 


(3.3) 

VA < 


G R, K(g + A z) > K(jf) 
particular, we have, 

A 

E(g + 

A 
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which implies that the Gateaux derivative of /C(g ) 


(3.4) 


dlCjg) = dJ(T{g),g) 

dg g 


Lemma 3.1. The mapping g -*■ T(g), from W(r<;) into H\Cl), has a Gateaux 
derivative ( dT(g)/dg)-z in every direction z E W(T C ). Furthermore, ( dT(g)/dg)-z = 
V (z) is the solution of 


(3.5) a{V,S) + c(u,V,S) n ,-{v,S) = 0 VS E H\Sl), 


(3.6) (V, R)r c = (z, R)r c Vi? E H~ l '\T e ). 

Proof. It is immediate from the linearity of (2.1)-(2.2). □ 

Now, we derive an optimality system from the first-order necessary condition 
(3.4). For each fixed g, the derivative dlC(g)/dg ■ z for every direction 2 E W(T C ) 
may be easily computed 

(3.7) 

■ z ss k 2 S f (Vg ■Vz + gz)dT+- f ( T - T d )V dT 
dg J v c 7 Jr. 

= k 2 £(V 5 , Vz)r c + K 2 S{g, z) Vc + ±(T - T d , V) a Vz E W(T e ), 

where for each 2 E W(r e ), V E is the solution of (3.5)-(3.6). 

Let (T,t) E H\n) x H~ 1/2 {T e ) be the solution of (2.1>(2.2) and let ($,r) € 
H 1 ^) x be defined as the solution of the adjoint problem 

(3.8) a(Z,*) + c(u,Z,$)n a + (r,Z)r e = ^{Z,T -T d ) v . VZ E H\ ft), 

(3.9) (W, $)r c = 0 VW E 

Setting 5 = $ in (3.5)-(3.6) and Z = V in (3.8)-(3.9), we have that 

(3.10) (t,z) r e = -(V,T-T d )r,- 

7 

Thus, from the necessary condition (3.4), we see that the optimal value of the control 
g satisfies 

(3.11) (V 5 ,V2)r c + (g,z) rc = -^( r > 2 0rc V 2 e ^(r c ). 

Collecting the above results, we obtain the optimal system 
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(3.12) a(T,S) + c(u,T,S) n ,-(t,S) re = (Q,S) n VS € H\n), 

(3.13) (T, R ) Fc - (g, R) Fe = 0 Vi? € 5" 1/2 ( r e ), 

(3.14) a(Z, $) + c(u, Z, $) n , + (r, Z) Pc = -( Z , T - r„)r, VZ € J5T l (12), 

7 

(3.15) (W, $) Fc = 0 We H~ 1/2 (T e ) 
and 

(3.16) (V ? ,Vz) re + ( 5 , 2 ) rc = -i-(r, ? ) re Vz£W( T e ). 

Integration by parts may be used to show that the system (3.12)-(3.16) constitutes 
a weak formulation of the boundary value problem 

(3.17) -Ki&T = Qi in fii, 

(3.18) -k 2 AT + (u • V)T = Qi + 2^(Vu + Vu T ) : (Vu 4- Vu T ) in fl 2 , 

(3.19) T = g on r e , 

(3.20) -«iA$ = 0 in fl 1? 

(3.21) — k 2 A$ - (u • V)$ = 0 in fl 2 , 
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4. Finite Element Approximation and Error Estimates 

A finite element discretization of the optimality system (3.12)-(3.16) is defined 
as follows. One may choose families of finite dimensional subspaces Vj* C H 1 ( fli), 
V£ C B’ 1 (Q 2 ) such that V r / l |r„ = These families are parameterized by 

the parameter h that tends to zero; commonly, this parameter is chosen to be some 
measure of the grid size in a subdivision of f2 into finite elements. Let V k = V”/ 1 U V k 
and 0 k = V fc |r e , i.e., O k consists of the restriction, to the boundary r e , of the 
functions, belonging to V h . For all choices of conforming finite element spaces, we 
then have that V k C H and 0 k C Next, let N h = F A |r e . Again, 

for all choices of conforming finite element spaces V h we have that N h c H l (T c ). 
Let Nq = N h fl VU(r e ). For the subspaces V k , V k , O h and N k , we assume the 
approximation properties: there exist an integer k and a constant C, independent 
of h, 7\, T 2 , t and g, such that 


(4.1) 

inf \\T l -T?\\ 1 <Ch m \\T 1 \\ m+ i,n l 

T* € V* 

VTi € H m+1 (ft!), 1 < m < k, 

(4.2) 

inf JlTt-TZlU <Ch m \\T 2 \\ m+ltC i, 

Tfevf 

VT 2 6 ff m+1 (fi 2 ), l<m<k, 

(4.3) 

mf ||t-t A ||_ 1/2 ,r c <C/i m ||t|| m _ 1/2 ,r c 

t * € O* 

Vt £ r c ), 1 < m < k 

and 



(4.4) 

inf || g - /|Ur c < C7i m - J+1/ %IU + i/ 2 .r c 


V<i £ VF(T C ), l<m<Jt, 0<s<l. 

A finite element algorithm for determining approximations of the solution of the 
optimality system (3.12)-(3.16) is as follows: seek T h £ V h , t h € 0 k , g h € N k , 
£ V h and r k £ O k such that 

(4.5) a(T k , S h ) + c(u, T h , S A )n, - ( t h , S A )r c = (Q, S h ) VS A £ V h , 

(4.6) (T\ R h ) Vc - (g k , R k ) Vc = 0 R k £ 0\ 

(4.7) K 2 6(g h ,K h ) rc + K 2 6(V s g k , V,K k )r c = -(K k ,r k ) rc VK h £ N k , 

(4.8) a(Z A ,* A ) + c(u,Z\$ A )n a + (Z k ,r k ) = ~(Z k ,T h - T d ) Vr VZ h £ V k 

7 
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and 

(4.9) (W\ $ A )r e = 0 W‘6 O h (r e ). 

The major task in this section is to obtain error estimates for the finite element 
solutions. It turns out to be convenient to apply the Brezzi-Rappaz-Raviart theory, 
even though our problem is linear. We introduce some spaces and operators, and 
verify the requirements of that theory. In the following discussion, the constants « lt 
« 2 and will be held fixed. Thus, the system (3.12)-(3.16) and (4.5)-(4.9) depend 
on the single parameter 7. 

Let X = JT(fl) x H- l ' 2 (T e ) x W(T e ) x E l (Q) x Y = x 

H~ 1/2 (T c ) x and Z = L 2 (fi) x L 2 (r e ) x L 2 (r ff ). Let the operator B E 

£(Y, X) be defined a s follows: B(Q,Q,P) = for ( Q,Q,P ) € Y and 

(T, t, g, $, f) € X if and only if 

(4.10) a(f,S) + c(u,f,S) a3 -(i,S)r e = (Q,S) n V5 6l‘(fi), 

(4.11) ( T,R) rc -{Q,R)r c = 0 VR e H~ 1/2 (T C ), 

(4.12) n 2 6(g, K)r c + '^sX)r ts = -{K,f) Tc VAT E W(r <: ), 

(4.13) a(Z, $) + c(u, Z, $)n, + (Z, f)r e = (^, P)n VZ £ 
and 

(4.14) (W,$)r e = 0 VW 6 J?~ 1/2 (n). 

Note that this system is weakly coupled. First, one may separately solve the prob- 
lems (4.10)-(4.11) for T and * and (4.13)-(4.14) for $ and f; then, one may solve 
the Laplacian problem (4.12) for g. 

Analogously, the operator B h 6 C{Y\X) is defined as follows: B h (Q,Q,P ) = 
(T h ,t h ,g h ,$ h ,f h ) for (Q, 0, P) € Y and ( f h ,t h ,g h ,$,f h ) 6 X h if and only if 

(4.15) a(f h ,S h ) + c(u,f\S k ) n ,-(t h ,S h )r c = (Q,S h )n VS h 6 V h , 

(4.16) (T k ,R h ) r c - (0, R h )r c = 0 E O h , 

(4.17) K 2 6(g h ,K h ) rc + M(V,S\ V.ff*) r e = -(K h ,f h )r c VA' h E N k , 


(4.18) 


a(Z\ * k ) + c(u, Z\ + (Z\ f h )r c = (Z\ P h ) a VZ* E V 
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and 

(4.19) (W fc ,i*)r e = 0 'iW h £O h . 

The system (4.15)-(4.19) is weakly coupled in the same sense as the system (4.10)- 
(4.14). 

Theorem 4.1. The second order elliptic problem (4.10)-(4.14) has a unique solu- 
tion belonging to X. Assume that (4.1)-(4.4) hold. Then the discrete second order el- 
liptic problem (4.15)-(4.19) has a unique solution belonging to X h . Let (T, i , g, $, f ) 
and (T h ,t li ,g h ,$ h ,f >l ) denote the solutions of (4.10)-(4.14) and (4.15)-(4.19), re- 
spectively. Then we have that 

(4.20) || f - f A || : + ||t - i*||_ 1/2 ,r e + II 9 - g h \\i,r c + 11$ - $*l|i 

+ ||f - f A ||_ 1/2i r c — 0 as h *• 0 

In addition, if(f,t,g,$,r) € X H lf2 {T c ) x H\T e ) x JJ 2 (fl) x H l ' 2 (T e ), then 

II T - T-^IU + ||f - f A ||- 1/2 ,r c + || g - g%,r c 

(4.21) + ||$ - $ A ||i + || f- •f A ||_ 1/2 ,r e 

< C'/i(||T , ||^r 3 ( n ) + ||$||/f3(n))- 

Proof. First, it follows from Lemma 2.1 that the two second order elliptic problems 
(4.10)-(4.11) and (4.13)-(4.14) each have a unique solution (T,f) and ($,r) belong- 
ing to 5" : (fi) x H~ 1/2 (T C ), respectively. From the Babuska’s theory, the discrete 
second order elliptic problems (4.15)-(4.16) and (4.18)-(4.19) each have a unique 
solution ( T h ,t h ) and ($ A , f h ) belonging to V h x O h , respectively. Moreover, we 
have that 

(4.22) \\T - r*|| x + ||t — ?||-i /3 . r . — 0 
and 

(4.23) ll$-$ A ||i + l|r-f A ||_ 1/2 . rc -0 

as h - 0, and if in addition (f,f) € H 2 {fl) x # 1/2 (r c ) and (*,f) E H 2 (Q) x 
H l ' 2 (T c ), we have that 

(4.24) ||T - T h ||i + ||f- PH-i/w. < CA||r||*, (n)> 

11$ - $ A ||i + ||r - r A ||_ 1/2 ,r, < C/i||$||*, (n) . 


(4.25) 
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Note that the problem (4.12) is a well known equation. Thus, we have that the 
problems (4.12) and (4.17) both have unique solutions, that 

(4.26) IliHkrc < C'l 1^1 |i/2,r c < C||$l|j»(n)> 

(4.27) I ! g - p*||i,r« < C{\\g - ^l|i,r e + Ik - f A ||_ 1/2 . r J Vp A € N 0 \ 

Using (4.4), (4.22) and (4.23) we then have that 

(4.28) llff- <kl|i,r c — *■ 0 as h -*• 0, 
and using (4.4), (4.24) and (4.25), we conclude that 

(4.29) ||p — p A ||i,rc ^ 

Then, (4.22), (4.23) and (4.28) yield (4.20), and (4.24), (4.25) and (4.29) yield 
(4.21). □ 

Let A denote a compact subset of R+. We define the operator Q from A x 
X to Y as follows: Q(y,(T,t,g,$,r )) = ( Q,Q,P ) for every ( Q,Q,P ) 6 Y and 
(7 , (T,t,p, $,r)) 6 A X X if and only if 

(4.30) (Q, 5) n = -(Qi,S) ai - (<? 2 , S) n , VS e 

(4.31) (0, R) Te = (- 9 , R)v c VR e ff- 1/2 (r e ) 

and 

(4.32) (P,Z) n = --(Z,T-T d ) r , VZ € 

7 

where Q 2 — Q 2 + 2/x(Vu + Vu T ) : (Vu + Vu T ). The operator Q is obviously of 
class C°°. The derivative of Q with respective to ( T,t,g , $,r), which we denote by 
Gx(l , (T,t,p, $,r)), can be defined as follows: 

(4.33) Qx(l, (T, t, p, $, r)) = (0, -p, ~^T) 

for every (7, (T,t,g, $, r)) 6 A x X. Furthermore, £?(7,(T,t,p, $,r)) e £(X,Z). 
Since A is a compact interval in R+ and the constant /c 2 is fixed, we see that Q and 
it’s first and second Frechet derivatives and all locally bounded maps. 

It is easily seen that the optimality system (3.12)-(3.16) is equivalent to 

(4.34) (T,t,p,$,r) + ££(7,(T,t,p,$,r)) = 0 

and that the discrete optimality system (4.5)-(4.9) is equivalent to 

(4.35) (T\ t\g\ r A ) + B h Q{ 7, (T A , t\g\ $\ r h ) = 0. 
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Now, having verified the hypotheses of the Brezzi-Rappaz-Raviart theory, we may 
use that theory to deduce the estimate 

||T - r*|| 1 + ||f - f*||_i/2,r c + l|<7 ~ <7*l|i.r c 

(4.36) + ||$ — $ A ||i + ||r — r A ||_ l / 2 ,r e 

< Ch(\\T\\ 6 , {n) + ||$||tf3(n))- 

We also have, from the theorem 2.5, applied to (3.12)-(3.16), the regularity estimates 

(4.37) ||r||^ (n) + !|i||i/2,r e + | |^tli.r c + ll^ll^(n) + ||T||i/ 2 ,r c 

< C(IIQI|o + IITrflli/’.r,)- 

The combination of (4.36) and (4.37) results in the following error estimates. 

Theorem 4.2. Let (T, t, $, r) be the solution of (3.12)-(3.16) and let ( T h ,t k , r h ) 
be the solution of (4.5)-(4.9). Assume that T, $ € also assume that (4.1)- 

(4.2) hold. Then, 

II T - T k ||! 4- ||* - t^ [ |— i/ 2 .r c 4- ||tf ~ ? h ||i,r c 

(4.38) + P - + || t - r h ||_ 1 / 2 ,r c 

<cfi(||giio + ||r rf || 1/ 2. r J, 

where C is independent of h, T, $ . 

We note that higher oreder estimates are possible if T is smooth in each sub- 
domain fix and fi 2 - 


5. Numerical Algorithm 

Let us consider the gradient method for the following minimization problem: 
Find g £ W(r c ) such that IC(g) := J(T(g),g) is minimized where T{g) £ H l { Cl) is 
defined as solution of (2.1)-(2.2). 

The classical Simple Gradient Algorithm proceeds as follows: 

Given g^; 

(5.1) define g ( n+ ^ = g ^ 1 d/C(g — ) recurs j ve ] v 

k 2 o dg 

Recall from §3 that for each fixed g, the derivative dJC(g)/dg ■ z may be computed 

(5.2) • z = k 2 6(- A s g + g,z)+ -( T - T d , V) r . Vz € W(T e ), 

dg 7 

where for each z £ W(r c ), V £ J ff 1 (0) is the solution of 

(5.3) a{V , S) + c(u, V, 5) n , = 0 VS £ H^Q), 
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V = z on T,. 


<9* 1 

(-^^)r c = ~(V,T-T d )r r . 

an 7 


(5.4) 

From (3.13), we see that 

(5.5) 

Thus, (5.1) may be replaced by 
for n = 0, 1,2, . . ., 

set </ (n+1) = gW - -L(K 2< J(-A,5 (n) + ? (n) )- K 2 ^-| re ) 


(5.6) 


= A 5 (») + I^!l 
,9 + 6 dn lr ‘’ 


where is determined from through the relations 
(5.7) a(T<- n \S) + c{u,T( n \S) = (Q,S) VS £ 


(5.8) T< n) = <7 (n) on r e , 
and 

(5.9) a(Z, $<">) + c(u, Z, $(">) = -(Z, T (n) - T d ) r , VZ € ff£(fl), 

7 


(5.10) 


$ (n) = 0 on r e . 


Therefore, we have the gradient algorithm results in the following iteration: 
Choose g 

for n = 1, 2, 3, . . . , solve for and from 
o(T<">,5) + c(u,r< n >,S) = (i Q,S ) VS € #£(fl), 

r (n) = 5 (n) on r c , 


and 

(5.11) i 

a(Z, *<">) + c(u, Z, $<">) = -(Z,T<"> - r d ) r , VZ G H l D ( ft), 

7 

$ (n) = 0 on T e , 
then solve for <^ n+1 ) from 


The convergence of the algorithm (5.11) is a direct consequence of the following 
lemma. 
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Lemma 5.1. Let K be a real-valued functional on a Hilbert space X with norm 
|| • | |x and scalar product (*, *)x. Suppose that there exist two constants m and M 
such that 

i) /C has a local minimum at a point x is of class C 2 in an open ball B centered 
at x } 

ii) Vu G B, V(x,y)eXxX, AC"(u) ■ (x,y) < M||x|U||y|| x , 

iii) Vu 6 B, Vx G X, IC"(u) • (x, i) > m\\x\\ 2 x . 

Let R denote the Riesz map, i.e. < f,x >= (R f,x)x for all x G X and all f G X' . 
Choose x(°) G B and choose a sequence {/>„} such that 0 < p. < p n < P' < 2m]M 2 . 
Then, the sequence {x^} defined by 

(5.12) *<") = xC" -1 ) - PnBXV” -1 )) for n=l,2,..., 

converges to x. Furthermore, if B — X and x is a global minimum, then the gradient 
algorithm converges to x for any initial value x^ 0 - 1 . 

Proof, see, e.g., [8] □ 

Theorem 5.2. Let (T^ n \ & n \g( n ^) be the solution of (5.11) and ( T , $,</) the solu- 
tion of (3.12)-(3.16). Then, if j6 is sufficiently large, g^ —* g and thus, T^ -*■ T 
in and — *• $ in H^Cl) as n — ► oo. 

Proof. In (5.11), we have the fixed parameter p = 1 /(k 2 S). For each g G VF(r c ), 
the second Frechet-derivative IC“(g) • ( z , w) may be computed by 

(5.13) X"(g) ■ ( z,w ) = k 2 6(Xw, Vz) Vc + K 2 6(w,z)r c + 1(17, V) r ,, 
where U G H 1 (£l) and V G H 1 ^) are the solution of 

(5.14) a(U,S) + c(u,U,S) = 0 VS G H^Q), 

(5.15) U = w on r e , 
and of 

(5.16) a(V, S) + c(u, V, S) = 0 VS G 


(5.17) 


U = z 


on r e . 
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One can easily have that ||£f||i < C| |tu||i,i' c and ||F||i < C||z||i,rc> where the value 
of the constant C depends only on fi. Then, 


(5.18) 


and 


(5.19) 


K"{g)-(z,w) < MIMkr.Mli.r, + — IMI1J.IWI1J. 

7 

= («,< + £)|Mkr.ll*lkr. 

«"(«)• (»,»)- + |V | 3 dT 

7 Jr, 

> Ml kl I* r.- 


Setting M = K^S + C/'f and m = k 2 6, we have, if jS > C/((\/ 2- 1)k 2 )> 2 mjM 2 > 
p = 1 /(k 2 S). The other hypotheses of Lemma 5.1 are easily shown to be valid. 
Hence, from that lemma, we obtain that 


(5.20) 


—> g in W(T C ) as n — ► oo. 


The desired convergence results follow from the a priori estimate (2.4). □ 


Of course, the gradient algorithm (5.11) is applied to the discrete equations. 
Then, we have two contribution to the errors in the computational solution, the 
discretization error T — T h and the iteration error T h — T h( ' n \ In a practical point 
of view, it is difficult to calculate A ,g^ in the last equation of (5.11). By using 
(3.24), we can substitute (5.11) by the following iteration: 


Choose and 

for n = 1,2,3,..., solve for T (n) and from 
a(T< n \ 5) + c(u, T< n) , S) = (Q, S ) VS € JTi(fl), 

r (n) = p (n) on r e , 


(5.21) 


and 


o(Z, $ (n) ) + c(u, Z, $ (n) ) = -(Z, T (n) - T d ) r, VZ 6 H l D (Q), 

7 


$ (n) = 0 on r c , 
then solve for g ( - n+1 ^ from 

0 c» + .) I _££ ™ 

9 ~ 9 6 V dn lr - dn 


O' 
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6. Numerical experiments 

Test 1 : We consider that the domain ft is the unit square (0, 1) x (0, 1) C R 2 , 
sub-domain fti = (0,1) x (0.75,1) and sub-domain ft 2 = (0,1) x (0,0.75). Let 
IV = (0.075, 1) x {0.75} C T w = (0, 1) x {0.75} and T c = {0} x (0, 0.75) (see Figure 
1 without the bump on the bottom boundary). 

The finite element spaces V j* and V* are chosen to be piecewise quadratic ele- 
ments on a triangle mesh such that V* = V£ on T w . We use the mesh size h — 1/12 
for all computation. Of course, calculations with varying mesh sizes were performed. 
In this paper, not being interested in the convergence history with varying mesh 
sizes, we do not report them. 

Now, we consider the following problem 


(6.1) 

-AT- 6.0 

on ft t , 

(6.2) 

- 2AT + (u • V)T = 

: 0 on ft 2 , 

(6.3) 

T = 1 + g 

on r„ 

(6.4) 

dT n 
— = 0 on 

an 

3ft\r c , 


where the velocity u is the solution of the Navier-Stokes equations 


(6.5) - Au + (u • V)u + Vp = 0 in ft 2 , 
the incompressibility constraint 

(6.6) div u = 0 in ft 2 , 
and the boundary condition 


(6.7) 

u = h 

on r e , 

(6.8) 

u = 0 

on iv u r 4 , 

(6.9) 

5^ = 0 and 
an 

ti 2 = o on r< 


where u = (uj,u 2 ) and h = (1.5y — 2y 2 ,0). To get approximate solutions for the 
Navier-Stokes solutions, we use the Taylor-Hood finite element on the domain ft 2 . 
Actually, we have simple solutions u = (1.5y — 2y 2 ,0) of the above Navier-Stokes 
problems. 



22 


MAX D. GUNZBURGER AND HYUNG C. LEE 



||T - 1 . 2 || 0 .r, 


J(T,g) 

7 = 1 , 6 = 1 
7 = 0.01, 6 = 0.003 

2.88669 x 10 - 2 
1.93423 x 10 " 3 

15.63322 

78.80300 

15.64765 

0.333120 


Table 1. The numerical results for Test 1. 


Note that since all data in (6.1)-(6.9) are sufficiently smooth in each domain fl, 
for i = 1, 2, we may assume that T G #*((2), s > 3 from the theorem 2.4. Thus we 
may use piecewise quadratic elements for the finite element approximation for the 
fast convergence with respect to ft. 

When g = 0 in (6.3), we say that the problem (6.1)-(6.9) is an uncontrolled 
problem. The numerical solution of the uncontrolled problem is shown in Figure 2 
and Figure 3 in which one can see that the temperature is above 2.0 on (0.3, 1) x 
{0.75} and even higher in the domain (0.3, 1) X (0.75, 1). 

Now, we try to get the desired temperature distribution along One can choose 
any reasonable desired temperature T d oa but we choose the parameter T d = 1.2 
on r, , thus we have 

(6.10) J(T,g) = L f | T- 1.2| 2 dT + 6 f (| 5 | 2 + |V j5 | 2 ) dT. 

*7 Jr c 

For the various choices of the parameters 7 and 6 appearing in the functional (6.10), 
the computations were performed. In this paper, we report the numerical results 
for the cases 

( 1 ) 7 = 6 = 1 , 

(2) 7 = 0.01 and 6 = 0.003. 

The costs are shown in the Table 6. 

In Figure 4-7, we plot the surfaces of the temperature T and adjoint state $ for 
each case. If one chooses the relatively small 7 6 , then one can have the relatively 
small value of ||T - 1.2|| 0 ,r,- 

Further reinforcement of our conclusions can be obtained from Figure 8 and 9 in 
which are found contour plots of the temperature T and adjoint state $. 

In Figure 10, we plot the approximate optimal control g h on the boundary T c . In 
Figure 11, we compare the temperature distribution on r„ in the uncontrolled case 
with the optimal temperature distributions in the controlled cases. 


Test 2 : We solve the problem (6.1)-(6.9) with h = (1.5y— 2y 2 , 0) on the domain fl 
which has a bumped boundary (see Figure 1 and 12). We assume that all parameters 
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Figure 2. The temperature surface plot for the uncontrolled problem. 



x 


Figure 3. The temperature contour plot for the uncontrolled problem. 
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Figure 4. The surface plot for the temperature T (7 = 8 = 1). 


Ci 

* 



Figure 5. The surface plot for the adjoint state $ (7 = 6 = 1). 
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Figure 6. The surface plot for the temperature T (7 = 0.01, 6 = 0.003). 



Figure 7. The surface plot for the adjoint state $ (7 = 0.01, 6 = 0.003). 
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II T - 1.2| 0 ,r» 

IMk 

J(T, g) 

7 = 1, s = i 

7 = 0.01, 6 = 0.002 

3.77877 x 10" 2 
3.13342 x 10" 3 

18.65837 

117.14224 

18.67727 

0.390956 


Table 2. The numerical results for Test 2. 


and data are the same as in Test 1. To get the approximate solutions for the Navier- 
Stokes equations, we also use the Taylor-Hood finite element on the domain fi 2 - 
We report the numerical results for the cases 

(1) 7 = « = 1, 

(2) 7 = 0.01 and 6 = 0.002. 

The costs are shown in the Table 6. We get the almost same results as in Test 
1 except that we need a little more control g on T c . Thus, even though the fluid 
flow is moderately complicated, given any 6 > 0, we can have 7 and 6 such that 
||T — Xdllo.rv < 6 when is sufficiently small. 

In Figure 12, we plot the temperature contour for the uncontrolled problem. In 
Figure 13-14, we have the contour plots of the temperature T and adjoint state $ 
for each cases. Finally, Figure 15-16 display the approximate optimal control g h 
along r c and the temperature distributions on T*, respectively. 

Remarks : For the case 7 = S = 1 , it was found that 10 — 15 iterations were 
sufficient to get the optimal control g. Since v = 1 and maximum velocity is 1, 
the control g affects the temperature distribution on T a very weakly. For the case 
that 7 6 is small, for example 7 S < 0.1, our gradient method does not converge. 
Thus, we need to adjust the iteration step size. In such case, we need a significant 
number of iterations. Thus, one may look for an efficient iteration algorithm. But 
the good news is that the iteration algorithm requires only one LU factorization and 
the same number of back and forward substitution as the iteration number , i.e., a 
comparable number of floating point operations relative to that required for solving 
the full coupled system (4.5)-(4.9). Of course we assume that h is sufficiently small. 






Uncontrolled 
(5 = 1 ,7=1 
6 = 0 . 01 , 7 = 0.002 
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Figure 16. The temperature distributions on T w . 
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